Marine  Physical  Laboratory 


AASERT:  Measurements  of  Wave 
Breaking  and  Dissipation  over  the 
Contental  Shelf 


Principal  Investigator(s)  W.  Kendall  Melville 


Supported  by  the 

Office  of  Naval  Research 

Grant  Number:  N0001 4-97-1 -0644 


Final  Report 


June  2003 


BEST  AVAILABLE  COPY 


University  of  California ,  San  Diego 
Scripps  Institution  of  Oceanogaphy 


20040809  088 


REPORT  DOCUMENTATION  PAGE 


Form  Approved 
OMB  No.  0704-0188 


Public  reporing  burden  for  tWs  collection  of  information  is  estimated  to  average  1  hour  per  response,  including  the  time  for  reviewing  instructions,  searching  existing  data  sources 
gathering  and  maintaining  the  data  need  ed,  and  completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any  other  aspect  of 
this  collection  of  information,  Including  suggestions  for  reducing  this  burden,  to  Washington  Headquarters  Services,  Directorate  for  information  Operations  and  Reports,  1 215  Jefferson 
Davis  Highway,  Suite  1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project  (0704*0108),  Washington,  DC  20503. 

1.  Agency  Use  Only  (Leave  Blank).  I  2.  Report  Date.  I  3.  Report  Type  and  Dates  Covered. 

July  1 1 , 2003  Final  Report 


4.  Title  and  Subtitle. 

AASERT:  Measurements  of  Wave  Breaking  and  Dissipation  over 

the  Contental  Shelf 

_____ 

W.  Kendall  Melville 


7.  Performing  Monitoring  Agency  Names(s)  and  Address(es). 

University  of  California,  San  Diego 
Marine  Physical  Laboratory 
Scripps  Institution  of  Oceanography 
291  Rosecrans  Street 
San  Diego,  CA  92106 


9.  Sponsoring/Monitoring  Agency  Name(s)  and  Address(es). 
Office  of  Naval  Research 
Department  of  the  Navy 
800  North  Quincy  Street 
Arlington,  VA  2221 7-5660 
Atten:  Louis  Goodman,  ONR  322 

11.  Supplementary  Notes. 


5.  Funding  Numbers. 
N0001 4-97-1 -0644 


Project  No. 
Task  No. 


enorming  organization 
Report  Number. 


9.  ^onsorj^/l^onitoring  Agency 


12a.  Distribution/Availability  Statement. 

Approved  for  public  release;  distribution  is  unlimited. 


12b.  Distribution  Code. 


13.  Abstract  (Maximum  200  words). 

The  objectives  of  this  project  were  to  develop  and  use  remote  in  situ  techniques  to  measure  the 
kinematics  and  dynamics  of  breaking  waves  and  relate  them  to  the  dynamics  of  the  wavefield  and  the 
surface  mixed  layer.  Of  particular  interest  has  been  the  dissipation  of  surface  wave  energy  by 
breaking. 


14.  Subject  Terms. 

surface  wave  breaking,  dissipating  wave  energy 


15.  Number  of  Pages. 


16.  Price  Code. 


17.  Security  Classification 
of  Report. 

Unclassified 


FINAL  REPORT 


AASERT:  Measurements  of  Wave  Breaking  and  Dissipation  over  the 

Continental  Shelf 

W.  Kendall  Melville 
Scripps  Institution  of  Oceanography 
University  of  California,  San  Diego 
La  Jolla,  CA  92093-0213 

phone:  (858)  534-0478  fax:(858)534-7132  email:  kmelville@ucsd.edu 

Award  Number:  N000 14-97- 1-0644  (2/1/97  -  1/31/02) 
http://airsea.sio.ucsd.edu 

GRADUATE  STUDENT  SUPPORT: 

This  AASERT  grant  was  used  to  support  the  participation  and  graduate  study  of  the  following  students: 

EricLuft,  MS  2001 
Ronan  Gray,  MS  2001 
Genevieve  Lada,  MS  2002 

The  major  scientific  accomplishments  of  the  parent  grant  are  described  below. 

LONG-TERM  GOALS 

The  long-term  goals  of  this  research  have  been  to  quantify  the  incidence  of  surface  wave  breaking  and 
its  role  in  dissipating  wave  energy,  which  is  then  available  to  generate  currents  and  turbulence  in  the 
upper  ocean. 

OBJECTIVES 

The  objectives  of  this  project  were  to  develop  and  use  remote  and  in  situ  techniques  to  measure  the 
kinematics  and  dynamics  of  breaking  and  relate  them  to  the  dynamics  of  the  wave  field  and  the  surface 
mixed  layer.  Of  particular  interest  has  been  the  dissipation  of  surface  wave  energy  by  breaking. 

APPROACH 

The  approach  revolved  around  the  use  of  aerial  imaging  from  light  aircraft  in  the  coastal  zone  and 
quantitative  image  sequence  analysis  using  techniques  of  Particle  Imaging  Velocimetry  (PIV)  to 
measure  the  kinematics  of  whitecaps.  Kinematic  measurements  were  used  to  investigate  aspects  of 
surface  mixing  including  estimates  of  the  rate  of  mixing  of  the  surface  layer.  In  conjunction  with 
simple  scaling  arguments  to  relate  the  kinematics  of  breaking  to  the  dynamics,  and  laboratory 
experiments  on  breaking,  the  remote  data  were  used  to  investigate  wave  dissipation  and  the  momentum 
flux  from  waves  to  currents.  In  situ  measurements  of  waves  and  currents  to  compare  with  the  remote 
measurements  remains  a  topic  of  continuing  research. 


WORK  COMPLETED 


This  project  has  focused  on  the  Shoaling  Waves  Experiment  (SHOWEX)  conducted  off  Duck,  N. 
Carolina  in  the  fall  of  1999  with  other  supporting  measurements  off  Oahu,  Hawaii,  and  California. 
During  SHOWEX  our  imaging  system,  the  Modular  Aerial  Imaging  System  (MASS)  was  flown  on  the 
Long-EZ  by  Tim  Crawford  (deceased)  of  NOAA.  Using  the  data  acquired  from  that  experiment  we 
developed  algorithms  based  on  both  PIV  and  simpler  techniques  to  analyze  the  image  sequences  to 
measure  A(c)  dc,  the  mean  length  of  breaking  fronts  per  unit  area  of  ocean  traveling  with  speeds  in  the 
range  (c,  c+dc)  (Phillips,  1985).  This  statistic  forms  the  basis  of  a  hierarchy  of  moments  that  can  be 
used  to  investigate  the  kinematics  and  dynamics  of  breaking.  Figure  1  shows  examples  of  the  image 
data,  the  traditional  whitecap  coverage  and  an  example  of  the  PIV  processing  on  a  single  whitecap. 

Two  moorings  were  deployed  to  measure  near-surface  turbulence  (ADV  and  Dopbeam  sensors)  and 
temperature  fields  and  bubbles  entrained  by  breaking  waves.  Two  bottom  ADCP/pressure  packages 
were  deployed  to  measure  waves  and  currents  on  the  shelf.  One  mooring  suffered  a  failure  during  a 
hurricane  and  was  redeployed.  One  of  the  moorings  was  inexplicably  destroyed  by  the  Canadian  Coast 
Guard  vessel  Creed,  a  participant  in  the  experiment,  during  good  weather  with  unlimited  visibility. 

RESULTS 

MASS  data  from  SHOWEX  has  been  analyzed,  written  up  and  published  by  Nature  (Melville  & 
Matusov,  cover  article,  May  2,  2002).  We  have  shown  that  during  SHOWEX  A(c)dc  increases  like  the 
wind  speed  cubed,  Uio3,  and  decays  exponentially  with  the  breaker  speed,  c.  Since  the  breaking 
statistics  depend  on  the  directional  wave  spectrum,  and  since  there  is  no  universal  wave  spectrum  over 
the  full  range  of  breaking  waves,  we  do  not  expect  any  universal  shape  for  A(c)dc.  Simple  physical 
and  dimensional  arguments  provide  interpretations  of  the  various  moments  of  A(c).  For  example,  the 
first  moment  is  proportional  to  the  rate  of  mixing  of  the  surface  area  by  breaking,  the  third  is 
proportional  to  the  volume  mixed,  the  fourth  is  proportional  to  the  momentum  flux  from  waves  to 
currents  and  the  fifth  is  proportional  to  the  wave  dissipation.  These  moments  are  shown  in  Figure  2. 
Note  that  as  you  move  to  higher  moments  the  mode  of  the  distribution  moves  to  larger  values  of  c. 
Since  c  scales  approximately  with  the  phase  speed  of  the  breaking  waves,  this  implies  that  while  small- 
scale  waves  contribute  most  to  the  length  of  the  breaking  fronts,  longer  waves  are  responsible  for 
momentum  transfer  and  dissipation.  Since  the  vector  c  and  not  just  the  speed  are  measured,  we  can 
examine  the  directional  distribution  of  A(c)  and  its  moments  can  be  measured.  Examples  in  Figure  3 
show  that  for  the  SHOWEX  data,  A(c)  and  its  moments  are  approximately  symmetric  about  the  wind 
direction. 

The  in  situ  data  from  the  surface  moorings  shows  that  the  dissipation  of  turbulent  kinetic  energy  (TKE) 
measured  by  the  Dopbeam  correlates  with  the  modulation  of  the  surface  waves  or  wave  groups.  An 
example  of  this  data  is  shown  in  Figure  4.  Analysis  of  the  Dopbeam,  ADV  and  temperature  data 
continues. 


IMPACT/APPLICATIONS 


During  this  project  we  have  shown  that  methods  of  PIV  developed  in  the  laboratory  can  be  used  to 
analyze  airborne  imagery  to  measure  the  kinematics  and,  with  supporting  assumptions,  the  dynamics  of 
wave  breaking.  These  results  are  novel  and  are  the  consequence  of  the  convergence  of  technical 
developments  in  the  areas  of  imaging,  position/motion  detection  (DGPS  and  inertial  systems)  and 
imaging  software.  While  we  have  concentrated  on  visible  imagery,  the  same  approach  could  be  taken 
with  IR  and  hyperspectral  imagery  to  study  the  kinematics  of  different  fields  including  temperature. 

The  in  situ  data  shows  great  promise  for  making  acoustical  measurements  of  ocean  microstructure 
directly  in  the  wavenumber  domain  (Veron  &  Melville,  1999).  This  avoids  the  use  of  Taylor’s 
hypothesis  to  move  from  the  frequency  domain  to  wavenumber  domain;  an  approach  that  has  been 
fraught  with  problems  in  the  surface  wave  zone  which  is  dominated  by  the  orbital  motion  of  the  waves. 

TRANSITIONS 

The  imaging  system  has  been  transitioned  to  the  measurement  of  breaking  and  spray  generation  in 
hurricanes  with  support  from  ONR’s  CBLAST  program. 

RELATED  PROJECTS 

Related  projects  are  described  on  our  web  page  at  http://airsea.sio.ucsd.edu. 
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Supplement:  N0001 4-97- 1-0644 

1 .  Number  of  Full  time  equivalent  graduate  students  (FTEGS)  supported  by 
parent  during  the  12  month  prior  to  the  AASERT  award  date  of  6/1/97. 

There  were  no  students  supported  in  the  parent  grant  the  year  prior  to  the 
AASERT  grant.  Parent  grant  was  N0001 4-97- 1-0277  (fund  31135A) 

2.  Total  dollars  of  parent  grant  $844,479 

3.  Total  FTEGS  supported  by  parent  grant  during  the  AASERT  grant  period. 

FTEG  graduate  student  - 10.5  months  at  the  GSR  rate  were  paid  for  G.  Lada 
and  F.  Veron.  This  equals  .875  FTEG  FTEG  undergrads  -  there  were  2382  hrs 
charged  of  various  students.  If  you  consider  that  a  student  would  be  half  time 
this  would  be  equivalent  to  a  little  more  than  2  FTEG.  Those  supported  were 
W.  Collins,  B.  Nahas,  R.  Gray,  B.  Lui,  D.  Norris,  N.  Serrano,  D.  West-Delgado 

4.  Total  AASERT  grant  FTEGS  and  undergraduate  students  supported  by  AASERT 
funds  during  the  grant  period. 

FTEG  graduate  student  -  there  were  3.29  FTEG.  G.  Lada  and  E.  Luft  were 
supported. 


Distance  off  shore,  m 


F  IG.Ia-c:  Images  showing  approximately  160  m  x  160m  of  the  sea  surface  for  mean  wind 
speeds,  Uio,  of  7.2,  9.8  and  13.6  m  s'1,  respectively,  in  SHOWEX.  Note  the  increasing 
density  of  whitecaps  with  Uio  and  their  random  shapes,  d  Whitecap  coverage  as  the  aircraft 
flew  offshore.  Note  the  large  fluctuations  at  spatial  scales  of  O  (1- 1 0)  km  that  may  be  due  to 
both  local  wind  and  wave  modulations,  e:  The  result  of  processing  a  single  whitecap  using 
the  PIV  technique  in  which  the  normal  velocity  of  the  boundary  of  the  whitecap  is  resolved. 
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FIG.  2a,  Binned  measurements  of  A(c)  weighted  by  Ui0'3  showing  that  the  average  length  of 
breaking  fronts  in  (c,  c+dc)  per  unit  area  of  sea  surface  increases  like  Uio3  and  decays 
exponentially  with  c:  A(c)  =  3.3x1  O^.e'0'640. 

b.  The  weighted  first  moment  of  A(c)  which  corresponds  to  the  fractional  area  swept  out 
(turned  over)  by  breakers  in  the  speed  range  (c,c+dc)  per  unit  time. 

c.  The  weighted  fourth  moment  of  A(c),  c4A(c)w,  which  corresponds  to  the  momentum  flux 
from  waves  to  currents  due  to  whitecaps,  with  the  subscript  “w”  referring  to  the  component  in 
the  wind  direction. 

d.  The  weighted  fifth  moment  of  A(c),  which  corresponds  to  the  energy  lost  from  the  wave 
field  due  to  breaking:  “wave  dissipation”. 


Figure  3a.  Distributions  of  A(c),  with  respect  to  the  northerly  and  easterly  components  of  c 
(cn.Ce),  along  with  the  average  wind  direction  (arrow).  Note  that  the  dominant  orientation  is 
close  to  the  wind  direction. 

b.  The  corresponding  momentum  flux  from  waves  to  currents  due  to  breaking  for  the  three 
wind-speeds.  Note  that  as  the  wind  increases  the  region  of  significant  momentum  flux 
becomes  approximately  symmetrical  about  the  downwind  direction  (c.f.  Fig.  2b). 

c.  The  corresponding  distribution  of  the  dissipation  with  (cn,Ce). 


FIG.4.  Top:  Inertial  estimates  of  dissipation  from  the  Dopbeam  acoustic  Doppler  near  the 
surface  in  SHOWEX. 

Bottom:  Simultaneous  measurements  of  the  vertical  component  of  orbital  velocity  measured 
with  an  ADV  on  the  same  mooring.  Note  the  corelation  between  the  modulation  of  the  wave 
field  and  the  dissipation,  especially  around  45-50  s. 
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